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ABSTRACT 

We present an analysis of BeppoSAX observations of three clusters of galaxies which are 
amongst the most luminous in the Universe: RXJ1347-1 145, Zwicky 3146 and Abell 2390. 
Using data from both the Low Energy (LECS) and Medium Energy (MECS) Concentrator 
Spectrometers, and a joint analysis with the Phoswich Detection System (PDS) data above 
10 keV, we constrain, with a relative uncertainty of between 7 and 42 per cent (90 per cent 
confidence level), the mean gas temperature in the three clusters. These measurements are 
checked against any possible non-thermal contribution to the plasma emission and are shown 
to be robust. 

We confirm that RXJ1347-1 145 has a gas temperature that lies in the range between 13.2 
and 22.3 keV at the 90 per cent confidence level, and is larger than 12.1 keV at 3cr level. The 
existence of such a hot galaxy cluster at redshift of about 0.45 implies an upper limit on the 
mean mass density in the Universe, f2 m , of 0.5. 

Combining the BeppoSAX estimates for gas temperature and luminosity of the three clus- 
ters presented in this work with ASCA measurements available in the literature, we obtain a 
slope of 2.7 in the L — T relation once the physical properties are corrected from the contam- 
ination from the central cooling flows. 
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1 INTRODUCTION 

The X-ray emitting gas in clusters of galaxies cools by the emis- 
sion of energy on a time scale which depends on the temperature 
and density of the intracluster medium (ICM). Where the density is 
higher, in the cores of the clusters, the cooling is more efficient and 
the temperature falls. To support the outer layers of gas, a subsonic 
flow of gas occurs towards the central region, producing a cooling 
flow that appears in X-rays as an enhancement of the central peak 
in the surface brightness profile and with a spectrum that presents 
multiple temperature components. 

The stronger the cooling flow, the higher the X-ray luminosity 
which arises from both the thermal energy lost by the cooling gas 
and the gravitational work done to maintain the gas at constant pres- 
sure. Therefore, to compare global cluster gas properties with theo- 
retical and numerical predictions, that, at present, cannot model the 
radiative cooling in detail, we need to account fully for the effects 
of cooling flows on the temperature and luminosity measurements 
(Fabian et al. 1994, Allen & Fabian 1998, Markevitch 1998). 

In this paper, we present data collected with the instruments 
onboard of the Italian-Dutch satellite BeppoSAX for three highly 
luminous (L bo i > 10 45 erg s" 1 ), hot (T gas > 10 keV), distant 
iz > 0.2) clusters of galaxies with large (>lOOOAf0 yr _1 ) 
mass deposition rates in their cores (Allen 2000): RXJ1347-1 145, 
Zwicky 3146 (Z3146) and Abell 2390 (A2390). Using data from 
both the Low Energy (LECS) and Medium Energy (MECS) Con- 



cetrator Spectrometer and a joint analysis with the Phoswich De- 
tection System (PDS) data above 10 keV, we show that it is possi- 
ble to constrain at a high level of significance the hard tail of the 
bremsstrahlung emission. 

All quantities quoted in this paper refer to the following cos- 
mological parameters: Ho = 50/i^q 1 km s~ 4 Mpc -1 , fi m = 1, 
Q A = 0. 



2 THE SAMPLE 

We describe here the characteristics of the three galaxy clusters 
analyzed in the present study. 

RXJ1347-1145 is a lensing cluster found in the ROSAT All 
Sky Survey (Schindler et al. 1995, 1997), with a slightly elongated 
X-ray structure around a very peaked central emission as it appears 
in the ROSAT HRI image. Analyses of the ASCA dataset provide 
temperature estimates ranging from 9.3^ { keV (Schindler et al. 
1997), when a single temperature model is adopted, to 26.4^2 8 3 
keV (Allen 2000), when a multi-phase gas is considered. 

Z3146 is the fourth ranked luminous cluster in the ROSAT 
band as reported in the Brightest Cluster Sample (BCS, Ebeling et 
al. 1998). Its massive cooling flow has been inferred from the strong 
optical lines and blue continuum exhibited by its central galaxy and 
has been studied with ROSAT HRI (Edge et al. 1994). Allen (2000) 
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quotes an ASCA estimate of the gas temperature corrected from the 
cooling flow of 11.3±i;f keV. 

A2390 is the sixth ranked luminous cluster in the ^OSAr band 
in the northern sky (BCS, Ebeling et al. 1998), has been widely 
studied in the X-ray waveband (Ulmer et al. 1986, Bohringer et al. 
1999) and is optically rich (Abell, Corwin & Olowin 1989, Yee et 
al. 1996). It presents evidence of strong and weak lensing (Pello 
et al. 1991, Pierre et al. 1996, Squires et al. 1996). For A2390, 
Bohringer et al. (1999) measure fcT gas ~ 9 keV when an isother- 
mal model is used to fit the ASCA data. Once they add a cooling 
flow component and consider a ROSAT PSPC-ASC4 joint fit, the 
temperature rises to ll.lt \% keV. Using the ASCA dataset only, 
Allen (2000) measures fcT gas = 10.13togg keV with a single 
isothermal model, that rises to W.S^^j 5 keV with a cooling flow 
component included in the spectral fit. 

In Table 1, we quote for these three clusters the exposure and 
net count rates measured in each instrument on board BeppoSAX , 
the Galactic absorption estimated in direction of the X-ray centroid 
determined from the respective ROSAT HRI image (see Fig. [j]) and 
their redshift from literature. Given the assumed cosmology, one 
arcmin corresponds to 0.408, 0.324, 0.277 hj Mpc, at the redshift 
of RXJ1347-1145, Z3146, A2390, respectively. 



3 SPECTRAL ANALYSIS 

We present the data from the three instruments onboard the Italian- 
Dutch satellite BeppoSAX (Boella et al. 1997): LECS (0.1-4 keV), 
MECS (1.8-10.5 keV) and PDS (15-200 keV). In Table 1, we 
present the details of the observations. 

Using XSELECT v2.0, we extract the spectra after screening 
for events without an attitude solution (generally 1 count for every 
1,000 good events) from the standard cleaned and linearized event 
file of LECS and the merged event file collected from units 2 and 3 
of the MECS. 

For all three clusters, we have selected circular regions with 
radii of 8 arcmin around the source in the LECS image and 6 arcmin 
in the MECS image. These regions should contain more than 95 
per cent of the flux observed in the MECS and about 90 per cent of 
that collected from the LECS, given the standard calibration files 
for point sources (see below for the calibration files adopted in the 
present analysis). 

Using similar regions in detector coordinates, background 
spectra were extracted from the MECS23_bkg . evt file, a 475 
ksec exposure of blank fields at high Galactic latitude. Following 
the procedure described in Fiore, Guainazzi & Grandi (1999), we 
have checked that the count rates in regions of the source field far 
away from the two calibration sources, the MECS strongback and 
any evident serendipitous source are consistent with the count rates 
measured in the same regions (in detector coordinates, DETX and 
DETY) in the blank-field. 

For the spectral analysis, we have made use of the cali- 
bration files released on September 1997 from the SAX-SDC, 
specifically the redistribution matrix file lecs_sep97 . rmf 
and mecsl_sep97 . rmf , and the ancillary response file 
lecs.130_124_8.sep97.arf and mecs23_6_sep97 . arf , 
for the LECS and MECS, respectively. Moreover, we have used 
the SAXDAS routine lemat to obtain both the RMF and ARF files 
for the three datasets simulating the extend emission of the cluster 
with n ~ 100 point sources each with an estimated contribution 
to the X-ray flux determined by the opportunely rebinned image in 
RAW coordinates. 
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Figure 1. Raw ROSAT HRI images with pixel-width of 5 arcsec for the 
three clusters in our sample, with adaptively smoothed contours, obtained 
using a Gaussian kernel that encloses a number of counts equal to five times 
the local background, overlaid. 
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Table 1. Summary of the BeppoSAX observations of the clusters in our sample. The Galactic absorption Nn is in unit of 10 20 cm 2 (Dickey & Lockman 
1990). References for the redshift: [1] Schindler et al. 1995, [2] Allen et al. 1992, [3] Yee et al. 1996. 

Cluster RA, Dec (2000) z 7V H date LECS (0.1^1 keV) MECS (2-10.3 keV) PDS (13-60 keV) 

hms, deg Exp. (ksec) 10 — 2 cts s~ 1 Exp. (ksec) 10 _2 ctss _1 Exp. (ksec) 10~ 2 cts s _1 

RXJ1347-1145 13 47 31.0,-11 45 11 0.4510 [1] 4.8 2000 Jan 25 28.6 5.3±0.1 73.2 9.9±0.1 36.5 4.0±2.6 

Z3146 10 23 40.0, +04 11 11 0.2906 [2] 3.0 1999Dec07 13.1 6.2±0.2 38.4 9.3±0.2 17.9 8.2±3.6 

A2390 21 53 36.8, +17 41 44 0.2279 [3] 6.8 1999 May 28 33.0 9.0±0.2 76.1 16.5±0.2 35.4 2.1±2.4 



The counts were grouped to a minimum of 25 per bin in order 
that the x 2 statistic could be applied in the fitting analysis. 

In the process of fitting the spectra simultaneously, we have 
left free the intercalibration constant of the spectral normalization 
in the LECS with respect to that in the MECS, used as reference. 
This value varies between 0.71 and 0.73 when an effective area ap- 
propriate for a point-source is used in fitting the LECS data. This 
is in agreement with the present absolute calibration of the two in- 
struments. When the more appropriate ARE file created with lemat 
is used, the intercalibration constant is about 1 for all the best-fit 
results. 

The PDS background-corrected spectra were retrieved from 
the on-line SDC archive and analyzed using the redistribution ma- 
trix file pds_256_aug97_f ix jan98 . rmf . We group the PDS 
counts according to the grouping file provided by SDC. The PDS 
vs. MECS normalization factor was fixed to 0.86. No differences 
are observed between PDS spectra that do or do not include correc- 
tions due to applying a crystal-temperature-dependent Rise-Time 
threshold. 

We fit the spectra using the photoelectric absorption model 
that adopts the cross-sections calculated by Balucinska-Church and 
McCammon (1992), in front of a MEKAL emission model (Kaas- 
tra 1992, Liedhal et al. 1995) in XSPEC (version 11.0.1, Arnaud 
1996). The absorption is fixed to the Galactic value quoted in Ta- 
ble 1 obtained from the 0.675° x 0.675° pixel HI map of Dickey & 
Lockman (1990). The results of the fitting analysis are presented in 
Table 2 and shown in Figure ^. As comparison, we quote in Table 2 
the results for a thermal fit to the MECS data only. 

In the same Table, we also present the results obtained from 
fitting two further models that include an additional contribution to 
the single-phase plasma: (i) a power-law that can take into account 
any non-thermal contamination from an unresolved active galaxy; 
(ii) an intrinsically absorbed cooling flow component which can 
model the emission from the central multi-phase gas (see also dis- 
cussion for the ASCA data in Allen & Fabian 1998). The aim in 
fitting these models is to check the robustness of the gas temper- 
ature estimates against any additional contribution to the thermal 
emission. We note, however, that the LECS, which provides data 
in the soft X-ray part of the spectrum and is used to constrain the 
intrinsic absorption and the normalization of the cooling flow com- 
ponent, has an energy resolution (FWHM) of about 11 per cent 
at 3 keV, that is comparable to ASCA Gas scintillation Imaging 
Spectrometers (GIS) but is poorer by a factor of 4 than the Solid 
state Imaging Spectrometers (SIS), which also has an effective area 
larger by a factor of 10 than the LECS. Therefore, any constraint 
on the modeling of the low-energy spectrum obtained through the 
"1T+CF" model will be less tight of that provided ASCA data (see 
Allen 2000). In our analysis, we placed lower and upper limit on the 
intrinsic cluster absorption of 10 19 and 10 22 cm -2 respectively. We 
flag as unc the values which are unconstrained within these bands 
by our analysis. 



The quoted pseudo-bolometric luminosity is estimated over 
the energy range 0.01-100 keV for the "1T+CF" model. 

Generally, all of the spectral fits provide an acceptable reduced 
X 2 ■ We discuss in the next subsection the results obtained from both 
the "IT" and "1T+CF" model. Regarding the "lT+pow" model, 
we note that only A2390 shows a constrained contribution from a 
power-law component with a luminosity less than 10 45 erg s -1 (2- 
10 keV) and a photon index of 2.69l 2 '3* consistent with the typical 
value for active galaxies (e.g. Nandra et al. 1997). Both RXJ1347- 
1 145 and Z3 146 do not show reasonable values to model the power- 
law component. For reference, we quote in Table 2 the best-fit re- 
sults obtained for the three clusters in exam with the "lT+pow" 
model and fixing the photon index to 2 (Nandra et al. 1997). 



3.1 Comparison with ASCA and ROSAT results 

All of the clusters in our sample were observed with both the ASCA 
(Tanaka et al. 1994) and ROSAT (Trumper 1983) satellites. We dis- 
cuss below the previous estimates of the gas luminosity and tem- 
perature in comparison to the results presented in this work. 

• RXJ1347-1145. Schindler et al. (1997) estimate a luminos- 
ity of 76 x 10 44 erg s" 1 in the 2-10 keV ASCA /GIS band. 
This measurement is about 20 per cent lower than the value esti- 
mated from the ASCA /GIS3 for the same dataset by Allen (2000; 
^44,2-iokcV ~ 93.5, where hereafter L44 is the luminosity in unit 
of 10 44 erg s _1 ). From the fit of the MECS spectrum alone, we 
measure an intermediate value between the two ASCA estimates of 
83.5 x 10 44 ergs" 1 . 

Schindler et al. (1997) also report a luminosity of 7.3 and ~ 21 x 
10 45 erg s" 1 with ROSAT /ERI in the 0.1-2.4 keV and bolometric 
bands, respectively. Whereas the extrapolation of our IT best-fit 
provides only ~ 70 per cent of the ROSAT flux in the 0.1-2.4 keV 
band, the bolometric value is in good agreement with the result 
obtained from the 1T+CF model applied to the joint analysis of 
the 3 BeppoSAX instruments and about 7 per cent lower than the 
estimate in Allen (2000). 

More controversial is the case for the gas temperature. While the 
estimate in Allen (2000) obtained from fitting simultaneously the 
0.6-10 keV SIS and 1-10 keV GIS data is in agreement with the 
best-fit value of about 14 keV within the 90 per cent confidence 
interval, the value of 9.3^'J keV (90 per cent confidence level) 
quoted in Schindler et al. (1997) can be excluded from the Bep- 
poSAX data at 99.99 per cent confidence level. In fact, from our 
analysis, the 3<r (99.73 per cent) lower limit is given by a tempera- 
ture of 12.1 keV. 

We show in Fig. ^| the residuals produced from the BeppoSAX 
data when the model of Schindler et al. (1997) is used. In particular, 
it seems that most of the residuals appear at energies below ~ 1 
keV and above ~ 8 keV, in spectral region not well covered from 
the selected energy range in their ASCA analysis of 0.7-9 keV. 

• Z3146. The estimated intrinsic rest-frame luminosity from the 
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Figure 2. (Left) Data and best-fit folded MEKAL model absorbed by Galactic absorption. We plot the result of the joint fit analysis of the LECS, MECS and 
PDS data and the corresponding residuals in terms of a. (Right) Confidence contours (A% 2 = 2.30, 4.61 corresponding to a confidence level of 68.3 and 
90 per cent, respectively, for 2 degrees of freedom) relative to the best-fit on the two interesting parameters, the metal abundance and the gas temperature. 
The solid line encircles the accepted values given a thermal plus cooling flow component; the dashed line represents the best-fit results obtained with a single 
thermal component; the dot-dashed line is obtained using a thermal component plus an absorbed power law (see Table 2). 



5 



Table 2. Best-fit spectral parameters. The errors are at 90 per cent confidence limit. The absorption is fixed to the respective Galactic value and the 
redshift to the optical estimate (cf. Table 1). The model used in XSPEC is phabs (mekal) (IT), phabs (mekal + zphabs (powerlaw) ) (lT+pow), 
phabs (mekal + zphabs (cfmodel) ) (1T+CF). The column "Prob" gives the null hypothesis probability that the fit is acceptable. The column "M-pow" 
quotes the deposition rate for the 1T+CF model, the photon index of the power law for the lT+pow model. The column "Z/2-iOkcV — ^bol" quotes either 
the luminosity in the 2-10 keV band for the fit on only MECS data or the expected contribution in the same band from any power law component (given the 
nominal values of the fit and upper limit at 90 per cent confidence level for a fixed photon index of 2) or the bolometric value calculated from the 1T+CF 
model. Note that unc means unconstrained. 



data 



Model 



kT 

keV 



abundance 

Z/Z Q 



Xl (d.o.f.) Prob 



A7V H 
10 21 cm" 



M-pow 



^2-10keV 
10' 



ibol 
44 erg s _1 



only MECS 


IT 


14 25+ 1 - 79 





51 


LECS + MECS + PDS 


IT 


14 48+ 1 - 76 





52 




lT+pow 


1 2 fi4+ 2 - 58 

-2.65 





46 




1T+CF 


15.88±2;66 





55 



RXJ1347-1145 

+0.19 
-0.17 
+0.19 
-0.18 
+0.22 
-0.19 
+0.21 
-0.19 



0.989(131) 
0.993 (207) 
0.990 (205) 
0.999 (205) 



0.52 
0.52 
0.53 
0.49 



103.41 475 , 
unc 



2 (fixed) 
< 1999 



83.5 

(14.2, < 30.2) 
204.0 



only MECS 
LECS + MECS + PDS 



IT 

IT 
lT+pow 
1T+CF 



7.26 
7.60 
7.44 
7.75 



+0.90 
-0.75 
0.92 
-0.77 
0.95 
-0.85 
3.47 
0.86 



Z3146 

Q no + 0.13 

r, OO+0.13 

33+°- 13 
U - JO -0.12 
,+0.14 
-0.12 



0.33 1 



0.904 (97) 
0.992 (141) 
0.976 (139) 
1.003 (139) 



0.74 
0.51 
0.57 
0.47 



> 976.8 
unc 



2 (fixed) 
< 1267 



32.4 

(49.4, < 106.8) 
69.6 



only MECS 


IT 


9 76+ ' 76 
a - ' u -0.68 





30 


LECS + MECS + PDS 


IT 


10 17+ - 77 
iu - i '-0.69 





29 




lT+pow 


a yq + 1.23 
°- ' -1.08 





33 




1T+CF 


10 67+ - 91 





31 



A2390 

+0.08 
-0.08 
+0.09 
-0.08 
+0.05 
-0.09 
+0.08 
-0.08 



0.991 (142) 
1.074 (263) 
1.039 (261) 
1.058 (261) 



0.52 
0.20 
0.32 
0.25 



30.8 



-18.6 

-13.1 



2 (fixed) 

r, + 308 

-337 



522 1 



34.7 

(7.8, < 11.6) 
78.3 




channel energy (keV) 

Figure 3. The best-fit model in second line of Table 1 in Schindler et 
al. (1997) is here considered: kT = 9.3 keV, N H = 10 21 cm~ 3 , 
Z = 0.32Zq. Only the normalization is left free to vary. The rebinning 
is done for presentation purpose. 

RASS (0.1-2.4 keV, Ebeling et al. 1998) is 26.5 x 10 44 erg s _1 . 
Using the same spectral model of Ebeling and collaborators, we 
calculate a I/44,o.i-2.4k c v of 20.0. Extrapolating the results of our 
IT model, we measure 25.8 x 10 44 erg s _1 . Thus, there is evidence 
for a slight mismatch in the spectrum shape between ROSAT and 
BeppoSAX , despite a general agreement in the normalization. 

In the 2-10 keV band, the fit of the MECS spectrum provides 
a luminosity of 32.4 x 10 44 erg s~\ 13 per cent lower than the 
estimate from ASCA /GIS3 in Allen (2000). 

Regarding the plasma temperature, our estimates are consistent 
within the 90 per cent confidence level with the results in Allen 
(2000). 



• A2390. Bohringer et al. (1999) fit both ROSAT /PSPC and 
ASCA /GIS data. They measure a luminosity of 27 (0.1-2.4 keV) 
and ~ 35 (2-10 keV) 10 44 erg s _1 . Ebeling et al. (1998) quote 
an intrinsic luminosity of 21.4 x 10 44 erg s _1 from the ROSAT 
All Sky Survey. Extrapolating to the 0.1-2.4 keV band the fit on 
our dataset, we obtain L44 = 23.7, 24.2 using the spectral model 
in Ebeling et al. (1998) and our best-fit IT model, respectively, in 
good agreement with the ROSAT estimates. 

Allen (2000) quotes a 2-10 keV ASCA /GIS3 luminosity of 
41.0 x 10 44 erg s _1 , that is about 17 per cent higher than both the 
estimate of Bohringer et al. (1999) and our measure from MECS 
data only of 34.7 x 10 44 erg s" 1 . 

The gas temperature quoted in both Allen (2000) and Bohringer 
et al. (1999), and obtained from analysis of ASCA data, are consis- 
tent within the 90 per cent confidence limits with the values mea- 
sured from the joint LECS-MECS-PDS fit. 

In conclusion, our results point to a general 10 per cent sys- 
tematic disagreement among the ASCA , ROSAT and BeppoSAX 
flux calibrations. It is worth noting that Grandi and Guainazzi show 
in the 1999 report on the BeppoSAX - ASCA intercalibration us- 
ing 3C273 that GIS3 provides a normalization of the continuum 
higher by a factor between 9 and 15 per cent than the correspond- 
ing value estimated in the same energy band with MECS for quasi- 
simultaneous data acquired in 1996 and 1998. 



4 A COOLING-FLOW CORRECTED L-T RELATION 

The relations between the measured (e.g. gas temperature and lu- 
minosity) and derivative (e.g. total gravitating mass) quantities de- 
termined from X-ray observations is a basic step in constraining 
cosmological parameters through the cluster mass function. 

The observed values of the intracluster gas luminosity and 
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temperature are correlated due to the fact that the luminosity de- 
pends on the cluster baryon fraction and the temperature is a direct 
measure of the potential well of the dark halo. 

In fact, the luminosity, L, goes like A(T)n 2 as r 3 ~ 
A(T)T 3/2 (1 + z) 3/2 (fi m /fi m (») 1/2 , where we have considered 
that (i) the gas fraction ratio is constant with redshift, (ii) the dark 
matter density evolves like (1 + jz) 3 , and (iii) the equation of the 
hydrostatic equilibrium, A/ to t ~ rT, holds. 

Considering that bremsstrahlung emission is the dominant 
emission process in clusters with kT > 2 keV, the cooling function, 
A(T), is then proportional to T 1//2 and, consequently, the bolo- 
metric luminosity is proportional to the square of the temperature, 
L ~ T 2 . However, this is only valid for a sufficiently hot cluster 
where the ambient luminosity and temperature are well determined 
above any contamination, like the excess in cool emission present 
in the core of cooling flow clusters (e.g. Fabian et al. 1994, Allen 
& Fabian 1998, Markevitch 1998). 

The three clusters presented in this work are among the most 
luminous known and are affected by the presence of the massive 
central cooling flows. Therefore, it is worth considering where they 
lie in the L — T plot once we take into consideration the best-fit re- 
sults of the "1T+CF" model. We add these results to a list of ASCA 
measurements of bolometric luminosity and temperature obtained 
from literature and corrected for the contamination from cooling 
flows, either by including a cooling flow model in the spectral anal- 
ysis (like we have done in the present analysis; cf. also Allen & 
Fabian 1998) or excluding the central emission to compose the an- 
alyzed spectrum (e.g. Markevitch 1998). 

The final sample is plotted in Fig. ^ and contains 64 en- 
tries. Of these, 3 are in common between our sample and the one 
of Allen and collaborators (A2390, Z3146, RXJ1347-1 145), five 
(A478, A1795, A2029, A2142, A2319) are analyzed both in Allen 
& Fabian (1998) and in Markevitch (1998). 

The physical values of the clusters in the combined sample 
range between 2.8 x 10 44 (A2657) and 2 x 10 46 (RXJ1347-1 145) 
erg s~ 4 in bolometric luminosity (the median of the distribution is 
3.0 x 10 45 erg s _1 ), between 3.5 (A1736) and 15.9 (RXJ1347- 
1 145) keV in temperature (median: 8 keV), between 0.04 (A2657) 
and 0.451 (RXJ1347-1 145) in redshift (median: 0.085). 

To investigate the behavior of the L — T relation, we use the 
following linear fit 

logF = loga + fclogX, (1) 

where Y = L x /(10 44 /i~ 2 erg s" 1 ) and X = T gas /(10keV), and 
we use the downhill simplex method implemented in the amoeba 
routine to minimize the following \ 2 merit function (Press et al. 
1992): 

a \- (logy ~loga-felogX) 2 
X (a, b)i = 



N 



+ b 2 ( 



(2) 



for each of the different datasets i considered in our sample. In the 
equation above, ey = (Tl/Lx and ex = or/Tgas- We assume a 
10 per cent uncertainty on the luminosity measurements (i.e., ey = 
0.1) and or equal to the maximum of the 68.3 per cent confidence 
asymmetric errors on the temperature. 

In particular, combining several different datasets, we follow 
the indications in Lahav et al. (2000) for minimizing a "weighted" 
merit function 



c; = ^iv l in (x 2 M)*) 



(3) 



where Ni is the number of data present in the dataset i. This for- 



malism allows us to consider some effective Hyper-Parameters, 
hp — Ni/x 2 (a,b)i, that provide a diagnostic on the systematic 
effects (more relevant for low hp) of that specific dataset in the 
overall modeling. Considering this, the BeppoSAX data are better 
modelled than the other two datasets (hp ~ 3 vs. 0.2). 

The errors on the best-fit parameters are estimated from 500 
bootstrap replications of the fit. In Table 3, we show both the best-fit 
parameters a and b and the values that characterize their distribution 
after 500 bootstraps (i.e., the median and the 16 and 84 percentile 
that corresponds to the lower and upper 1 a level in a Gaussian 
distribution). 

Here we note that the computation of the physical properties 
presented in the present work follows the technique used in Allen 
(2000) and differs from the analysis done by Markevitch in the fact 
that we try to model spectrally the emission from the core, whereas 
Markevitch cuts the core emission from the spectral analysis and 
extrapolates the contribution from the core to the total luminosity 
from the external emissivity. To take into account the different way 
in which the bolometric luminosity is estimated, we leave free the 
normalization, a, in the L — T relation in eqn. [l] when we perform a 
linear fit in each dataset. Only the slope parameter is linked among 
the three datasets. 

Combining the three datasets, we observe a slope in the L — T 
relation of about 2.7. From the distribution of the best-fit values 
after bootstrap resampling, we are still ~ 4<r far (too steep) from 
the value of 2 predicted from the simple gravitational collapse of 
plasma in a cluster dark matter. 

This result is consistent, within la, among the three different 
datasets here considered and with the Bayesian analysis applied to 
the Allen (2000)'s and Markevitch (1998)'s sample by Reichart, 
Castander & Nichol (1999). 



5 DISCUSSION AND CONCLUSIONS 

From the BeppoSAX observation, we can constrain the ambient 
temperature in the clusters of galaxies RXJ1347-1 145, Z3146 and 
A2390 with a relative uncertainty of (17/41), (11/41), (7/9) per cent 
in the (lower/upper) boundary (at the 90 per cent confidence level) 
respectively, when a cooling flow model is included in the spectral 
analysis. 

Although the instrument characteristics of BeppoSAX do not 
allow us to constrain firmly any absorption intrinsic to the clus- 
ter atmosphere and/or the normalization of the emission from the 
central cooling gas, the measurements of the ambient gas tempera- 
ture are in agreement, within the 90 per cent confidence level, with 
the results from the analysis of ASCA data in Allen (2000). In par- 
ticular, the plasma temperature in RXJ 1347-1 145 appears well de- 
fined between 13.2 and 22.3 keV (90 per cent confidence level) and 
above 12.1 keV at the 3a level. This measurement of the temper- 
ature is in good agreement with an independent estimate of 16.2 
(± 3.8, lcr level) keV obtained from the observed X-ray gas den- 
sity and the measured Sunyaev-Zeldovich (1972) distorsion of the 
cosmic microwave background in the direction of RXJ 1347- 1145 
(Pointecouteau et al. 1999; see also Komatsu et al. 1999). 

Following the arguments of many authors (e.g. Bahcall & Fan 
1998, Donahue et al. 1999, Henry 2000) the existence of such a hot 
cluster at redshift 0.45 is highly unlikely in an Einstein - de Sitter 
universe. We can estimate this probability, calculating the expected 
evolution of collapsed structures with a given minimum mass ac- 
cording to the Press-Schechter (1974) formalism. Integrating the 
number of collapsed objects with gravitating mass above that cor- 
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Table 3. Best-fit results on the luminosity-temperature relation for the data in our sample with median values of L = 30.0 X 10 44 erg s -1 and T = 8.0 
keV. The first two columns state the sample analysed and the number of objects in exam. The third column indicates the number of objects considered in 
each subsample and the measured \ 2 for the global best-fit. The scatter in the fourth column is defined as ^(log V — loga — 6 log X) 2 /TV] 1 / 2 , where 
Y = L x /(10 44 /i5 2 erg s" 1 ) and X = T gas /(10keV). Note that the scatter along the X -axis can be estimated as a\ n x = ""in y /b.ln the last column, we 
quote the best-fit slope after applying the technique described in Section 4, and the median and lcr error after 500 bootstrap replications. 



cut 


N data 


N sample (\ 2 ) 


scatter 


b 


only BeppoSAX data 


3 


3 (0.9) 


0.327 


2.21 (1.50 ± 1.13) 


only Allen (2000) 


30 


30(141.3) 


0.512 


2.41 (2.41 ±0.35) 


only Markevitch (2000) 


31 


31 (108.7) 


0.240 


2.78 (2.78 ±0.16) 


all 


64 


3 (1.0), 30(146.8), 31 (109.9) 


0.437 


2.69 (2.76 ±0.17) 


L > 10 X 10 44 ergs" 1 


53 


3 (1.0), 30(141.5), 20(71.0) 


0.420 


2.45 (2.57 ±0.29) 


T > 6 keV 


46 


3 (1.2), 25 (86.0), 18(61.6) 


0.388 


2.72 (2.87 ±0.30) 




Figure 4. The observed luminosity-temperature relation from BeppoSAX 
and ASCA data with the best-fit result obtained for the 64 clusters in our 
sample. The solid dots represent the BeppoSAX results for RXJ1347-1 145, 
Z3146 and A2390. The squares are from Markevitch (1998). The triangles 
are from Allen (2000). All the estimated temperatures are corrected for the 
presence of the cooling flow. The error bars are 90 per cent confidence limit. 
The solid, dotted, dashed lines are the best-fits obtained fixing to a common 
value the slope and leaving free the normalizations for the objects in the 
BeppoSAX , Allen (2000) and Markevitch (1998) samples, respectively. 

responding to a temperature of 15.9 keV over the redshift range of 
0.4 — 0.5, and requiring at least the existence of RXJ1347-1 145, 
the single detection provides a stringent upper limit of O m < 0.43. 
Considering also the lower and upper end of the range of the ac- 
ceptable gas temperature of 13.2 and 22.3 keV, with the 90 per cent 
uncertainty on the single detection, we constrain the cosmological 
parameter fl m between 0.25 and 0.70 (cf. Fig. ^). 

Future Chandra and XMM observations of these highly lumi- 
nous clusters of galaxies at intermediate redshift will resolve spa- 



Figure 5. The expected number density of clusters with gas temperature 
larger than 13.2, 15.9, 22.3 keV (from right to left in the figure, respectively; 
cf. Table 2 for the constraint on the gas temperature of RX11347-1 145) in 
the redshift interval [0.4, 0.5]. The dashed lines represent the distribution in 
a "f2 m + = 1" Universe (thick) and in a open Universe (thin). The error 
bars, around the observed value of 1 cluster with the given characteristics, 
are at the 90 per cent level of confidence, according to the Poisson single- 
sided distribution tabulated in Gehrels (1986). 

tially the profiles of the gas temperature, abundance and deposition 
rate. 
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